Incompressibility of polydisperse random close packed colloidal particles 
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We use confocal microscopy to study a random close packed sample of colloidal particles. We intro- 
duce an algorithm to estimate the size of each particle. Taking into account their sizes, we compute 
the compressibility of the sample as a function of wave vector q, and find that this compressibility 
vanishes linearly as q ^- 0. The particle sizes must be considered to calculate the compressibility 
properly. These results also suggest that the experimental packing is hyperuniform. 

PACS numbers: 82.70.-y, 61.20.-p, 64.70.pv, 64.70.kj 



The random packing of objects has been studied sci- 
entifically for nearly a century [l|, [2[; see Ref. Q for 
a review. This problem is often termed "random close 
acking" (rep) or "maximally random jammed packing" 
^_^ . Important recent work has focused on the packing of 
highly polydisperse systems Q, ellipsoids Q, and tetra- 
hcdra [7| , but the simplest packing problem is the packing 
of monodisperse spheres. In the past decade, simulations 
studying monodisperse spheres have generated large rep 
configurations with 10^ — 10^ spheres @, Q- These simu- 
lations enable study of density fiuctuations at very large 
length scales, or equivalent ly, small wave vectors q. They 
find that the static structure factor S{q) approaches zero 
linearly as q ^ 0, that is, S{q) ^ q for small q. This find- 
ing has been termed "hyperuniformity" [8[. One corol- 
lary is that the sample is incompressible, as the isother- 
mal compressibility x in simple liquids can be found from 
pksTx = S(0) where p, fcs, and T are the mean density, 
Boltzmann constant, and temperature. These observa- 
tions of close-packed samples are in contrast, for exam- 
ple, with simple liquids for which 5(0) > [1^. The 
existence of hyperuniformity has been seen in a va riety 
of systems, sec for example discussions in Refs. [3, uA\ - 
In general, long wavelength density fluctuations are im- 
portant for diverse fields including critic al p henomena 
[13 j and the shear fiow of glassy materials [IjI • Likewise, 
understanding random close packed samples is relevant 
for understanding liquids, glasses, biological systems, and 
granular materials [l|, y, [Ij] . 

In 2010 we published an experimental study of a ran- 
dom close packed sample of colloidal particles, observed 
with confocal microscopy [la] . Our data set was the posi- 
tions of more than 500 000 slightly polydisperse particles 
[i&j . and we found that S{q ^ 0) > 0, implying that the 
experimental sample was compressible and not hyperuni- 
form. A 2010 simulation of a binary sample found similar 
results [l7[. These results seem to demonstrate random 
close packed samples that are not hyperuniform. How- 
ever, in 2011 two groups showed that in polydisperse sam- 
ples, careful consideration of the individual particle sizes 
recovers hyperuniformity and incompressibility |10l . Illl ] . 
In particular, Berthier et al. showed how to compute the 
isothermal compressibility when the individual particle 



sizes are known, and demonstrated that sainples with 
5(0) > nonetheless can be incompressible [l3|- They 
examined data from a two-dimensional granular experi- 
ment and confirmed that x(0) = 0. The reason 5(0) > 
in polydisperse systems is because density fluctuations 
are coupled to composition fluctuations, but such sam- 
ples can still be incompressible and hyperuniform. 

In this article, we describe a method to determine each 
particle size from microscopy observations of a random 
close packed sample of colloidal particles. We use nu- 
merically generated packings to confirm that our method 
accurately determines the particle radii. Analyzing our 
experimental data using the method of Berthier et al. 
[lO|, we confirm that our experimental system is hype- 
runiform and incompressible. We additionally note an 
anticorrelation between the local polydispersity and lo- 
cal ordering. 

As we use the analytical method introduced by 
Berthier et al [lO[, we briefiy summarize their method 
here. They consider a wave vector dependent isother- 
mal compressibility xil) which is related to the structure 
factor of a monodisperse sample by pksTxiq) ~ S{q). 
They then derive an exact formula relating xil) ^-i^d 
S{q) for a polydisperse sample, although the formula 
is "conceptually and computationally difficult" to eval- 
uate [lO|. Thus, they derive a series of approximate 
formulas, of which the first order approximation is suf- 
ficient for samples of low polydispersity such as ours. 
To start with, they define single-particle density fields 
Pi(q) = exp(iq • rj) where rj is the position of parti- 
cle i. They also define the size deviation of particle i as 
Ej = (oi — a)/a, where at is the radius of particle i and a is 
the mean radius. (Note that \/Jefj = p defines the poly- 
dispersity p of a sample.) These e^'s are the smaU param- 
eters used in the approximation. Using these variables, 
they define a 2 x 2 matrix S{q) with elements S'"^{q) = 
i{e"iq)e^-q)), with u,v e 0, 1, e«(<7) = J:l,e^p,{q), 
and N is the total number of particles. The matrix ele- 
ments can be used to provide a first order approximation 
Xi(g) as pkeTxiiq) = 5°" - [S°^]^/S^\ They confirm 
that xi(0) ~ in cases for which the sample polydis- 
persity is less than 10%, while 5(0) 7^ for those cases. 



Their results suggest that random close packed systems 
are hyperuniform and incompressible even when the sam- 
ple is polydisperse |ig |. 

In our prior work, we used colloidal particles to gen- 
erate a random close packed sample, and imaged this 
with confocal microscopy. We reprise the key experimen- 
tal points here; a more detailed experimental discussion 
is in Ref. [151 . We use sterically stabilized poly(methy 
methacrylate) (PMMA) particles [l^ with a = 1.265 fim. 
Previously we reported that these particles had a polydis- 
persity of ~ 5% [l^; below, we determine that the true 
polydispersity is 6.7%. The PMMA particles are sus- 
pended in a solvent mixture that is slightly lower density 
than the particles. The sample is mixed and then the par- 
ticles are allowed to sediment until they are close packed. 
We use a confocal microscope to take clear images deep 
inside our dense sample [13- Overlapping images are 
taken, with total volume 492 x 514 x 28 /^m^. Within 
this volume, particles are identified within 0.03 /ini in x 
and y, and within 0.05 /xm in z [l9l . l2l| . The total data 
set contains 543 136 particles jla |. 

The average particle size a is obtained from the posi- 
tion of the first peak of the pair correlation function [1^ . 
It is difficult to determine subtle size differences between 
individual particles from microscopy due to diffraction. 
However, obtaining the positions of each particle can be 
done fairly accurately. A large particle will be slightly 
farther from its neighbors as compared to a small parti- 
cle, and we use this idea as a starting point for an esti- 
mation method for each particle size. 

Given that our sample is jammed, each particle must 
be in contact with several of its neighbors. In fact, a nu- 
merical simulation of random close packed monodisperse 
particles showed that each particle contacts with at least 
6 particles [20|. When particle i contacts with particle 
j, the separation between these two particles is given by 
Tij = Oi + Oj, where a^ and aj are their radii. The average 
of rij over all neighbors j is given by {rij)j = ai + {aj)j. 
Next, consider separations rjk between particle z's con- 
tacting neighbors j and contacting neighbors k of those 
particles. Again, we take an average of rjk with respect 
to particles j and fc, giving {{rjk)k)j = {aj)j + {{ak)k) 



Then we subtract {{rjk)k)j from {rij)j, leading to 
= {{ak)k)j + {nj)j - {{rjk)k)j- 



(1) 



We choose the 5 nearest particles from particle i as 
the particles j, assumed to be in contact with particle i, 
and likewise for each particle j we identify its five clos- 
est neighbors for the particles k. For each particle j, 
one of its neighbors k should be particle i, leading to 
an overcounting in the average: {{ak)k)j ~ (l/5)ai + 
(4/5)((afe)fe^i)j. Likewise, {{rjk)k)j = (l/5)(rij)j + 
(4/5)((rjfc)fc^i)j from the same overcounting of particle 
i. Using these results, we obtain 



ai = {{ak)k^i)j + {nj)j - {{rjk)k^i)j 



(2) 



a 


'.•y' 


1.1 


.Jr 




.^^fJF 


V^-" 


jr 


0.9 


;r 




0.9 1.0 1.1 



0.00 



0.00 0.04 0.08 0.12 
P 



FIG. 1: (Color online) (a) Scatter plot of the calculated ra- 
dius al from our method (Eqn. [2| as a function of the given 
radius a'g using data from a simulated packing with polydis- 
persity 7%. The solid line corresponds to al = ag. (b) The 
particle size uncertainty Aa found by analyzing simulation 
data from packings with a given polydispersity, both without 
noise (circles) and with noise added to the particle positions 
(triangles). The dashed line corresponds An — p. 



accurate values for a^. In this way at is found from the 
mean particle size and the particle separations, which are 
obtained directly from microscopy. 

To validate our method, we simulate polydisperse rep 
samples using the algorithm of Refs. [23, 123 ■ We use 
512 particles with mean radius a = 1 and polydispersity 
from 0.01 to 0.12, generating 5 independent configura- 
tions for each polydispersity. The particle size distribu- 
tion is a Gaussian. Using the simulated position cen- 
ters, we calculate the radii of the particles a^ by our 
method. Figure [U^a) shows a scatter plot of a* as a 
function of the given radii a* from a simulation with 7% 
polydispersity. The calculated radii are located around 
a 



a* . We define the uncertainty of the size estimation 



To compute {{ak)k^i)j, we use a as an initial guess for 
the particle sizes, and then iterate five times to get more 



as Aa = i/([(a^ — al)/a^g]'^)i. Aa is plotted as a func- 
tion of polydispersity p as circles in Fig. [TJb) . We find 
Aa w p/6. The polydispersity of a^ matches that of at. 

One experimental complication is that there is an un- 
certainty in the position of each particle. In our experi- 
ment, the uncertainties are 0.024a in x and y and 0.0395a 
in z. We add this positional uncertainty to the true sim- 
ulated positions, and then redetermine the particle radii. 
As expected, this increases the uncertainty Aa of the final 
radii, shown by the triangles in Fig. [IJb). Aa increases 
by ^ 0.01 compared to the case without positional noise. 
Positional noise is fatal when the polydispersity is less 
than 0.02, but otherwise our method results in more ac- 
curate radii even in the presence of noise. 

Next, we estimate each particle size of our experimen- 
tal data with our method. Given that Eqn. [2] requires 
information about both a particle's nearest neighbors 
and also second nearest neighbors, only particles suffi- 
ciently far from the edges of our images have accurate 
sizes. We modify our algorithm slightly for the experi- 
mental data as follows. We find the coordination number 
Zi of each particle, the number of neighboring particles 
within a distance 2.8a (the first minimum of the pair cor- 
relation function) [l5| . From the particles in the interior 
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FIG. 2: (Color online) Probability of particle sizes in our 
experimental sample. The average size is 1.265 /im and poly- 
dispersity is 6.7%. 



of the sample, we find the average coordination number 
z w 12. Then, for every particle, we estimate the num- 
ber of touching neighbors T^ = bzi/ 12 where we round 
Ti to the nearest integer. For particles at the edge of 
the imaged volume, T^ < 5 as not all of the neighbors 
are imaged. Then for each particle, when averages over 
contacting neighbors j are done in Eqn. [2J these aver- 
ages are over the Ti nearest neighbors. After iterating 
Eqn.[2]to find all radii, the edge particles are removed by 
cropping the data to a volume of 440 x 461 x 14.2 iim^, 
containing 217 816 particles. 

Based on these particles with their calculated sizes, the 
volume fraction of this sample is found to be (/) = 0.647 
± 0.007, where the uncertainty of (f> is due to the uncer- 
tainty in determination of each particle size. Figure [5] 
shows a distribution of the estimated particle sizes. This 
sample has a polydispcrsity of 6.7%. Given this mea- 
sured polydispersity, Fig. [IJb) shows that Aa w 0.023 
(corresponding to aAa = 0.03 /xm). The experimental 
distribution is not a Gaussian and this is not an artifact 
of our method, as a simulated Gaussian size distribution 
with positional noise leads to a measured Gaussian size 
distribution. 

Using our estimated particle sizes, we now study the 
wave vector dependence of the compressibility xo('?) 
and Xi{l) of our experimental data. Figure [3] shows 
pkBTxoiq) and pkBTxi{q)- Our experimental data do 
not obey periodic boundary conditions, and the effect of 
the boundaries appears near q = Q. xo{l) and Xi{l) a-^^ 
independent of the choice of Fourier window functions 
for qa/i: > 0.2. Thus we do a linear fit to pkBTxoiq) 
and pkBTxiiq) in the region 0.2 < qa/ir < 0.5, shown 
as the lines in Fig. |3l both functions have linear behav- 
ior in this region. We find pkBTxi{0) = 0.002 ± 0.004, 
while pkBTxo{0) = 0.049 ± 0.008 as reported previously 
[15| . The uncertainties are due to the uncertainties of 
particle positions and sizes, and the choice of the fitting 
range. Our observation that Xi(<?) ~ 9 shows that long 
wavelength density fluctuations are suppressed. This is 
consistent with the observations of Berthier et al. and 



FIG. 3: (Color online) pkBTxo{q) (no approximation) and 
pknTxiiq) (first order approximation of Ref. [ly|), from the 
experimental data. Square symbols correspond to pkBTxo{<l), 
which is proportional to S{q) at small q. Circle symbols cor- 
respond to pkBTxi{q)- The lines are linear fits to the data 
for 0.2 < qa/Tv < 0.5. 
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4; (Color online) The number of ordered neighbors No 
as a function of the local polydispersity pL in the experi- 
ment (open circles) and simulations (bulk polydispersities as 
indicated in legend). The line is the average of A^o from 75 
simulated systems with bulk polydispersity from 1% to 12%. 



show that our system is incompressible and likely hype- 
runiform [l0|. 

Our data let us consider a new question, the relation- 
ship between local environment and local ordering. It 
is known that crystallization occurs in samples with low 
polydispersity {p < 0.08) [2J-[2g. However, crystal nu- 
cleation is a microscopic phenomenon, that is, the crystal 
nuclei do not necessarily "know" the bulk polydispersity. 
We can use our data to investigate the relationship be- 
tween local ordering and local polydispersity. 

We define the local polydispersity p\ as 



Pl 



\/((a„ - ai)'^)/ai 



(3) 



where (a„) is the mean radius of the nearest neighbor par- 
ticles of particle i. The nearest neighbors of a particle are 
defined as those with centers separated by less than 2.8a 
p^ . We calculate the bond order parameter d™ to quan- 
tify how the local structure compares between neighbors 
i and n [27H29[. Two neighboring particles are termed 



"ordered neighbors" if dg" exceeds a threshold value of 



0.5 |27H29| . We then count the number of ordered neigh- 
bors Nl around each particle i. N^ ^ corresponds to 
random structure around particle i, while N^ > 7 means 
that particle i is in a crystalline environment [29[. 

Figure |4] shows that the local polydispersity p\^ has a 
strong influence on local order N^. The open circles show 
the result from our experiment. Particles with low p^ arc 
more ordered than particles with high p^: that is, there is 
a tendency for particles to order when the central particle 
size ai is similar to its surrounding neighbors. A similar 
result is found from our simulated packings (closed sym- 
bols and solid line in Fig.|4|), where the local polydisper- 
sity predicts local order independent of the global poly- 
dispersity. The agreement between the simulations and 
the experiment is striking, especially given that the sim- 
ulation corresponds to an extremely fast quench, whereas 
the experimental quench allows time for particles to re- 
arrange [ia|. Note that these conclusions are unchanged 
when Ui in Eqn.[3]is replaced by (a„), although the trend 
shown in Fig.|3]is less pronounced. Our results are consis- 
tent with the prior knowl edg e that polydispersity affects 
the ability to crystallize |24| - [2a |. but this is the first ex- 
amination we are aware of showing how polydispersity 
can have a local influence on crystallization. It suggests 
an intuitively reasonable idea, that in a moderately poly- 
disperse sample, crystalline nuclei are more likely to form 
from locally monodisperse patches. 



We note that the observed polydispersity of our sam- 
ple (6.7%) helps explain a discrepancy we noted between 
our observations [la] and those of DuUens et ai, who also 
studied dense suspensions of sedimenting particles with 
similar sedimentation rates [30[ . They observed that par- 
ticles formed crystals in all cases [30] , while our particles 
pack randomly. Their samples had a polydispersity of 
5%, while our sample is 6.7%. Crystal nuclcation is sen- 
sitive to polydispersity in this range [2g] and this likely 
explains why our sample avoids crystallization, and why 
the samples of Dullens et ai crystallized. 

To summarize, we have presented a method to esti- 
mate the sizes of individual colloidal particles from ex- 
perimental knowledge of only their positions, and relying 
on the fact that the sample is close-packed. Numerical 
simulations confirm that our method is robust even in 
the presence of realistic experimental noise. Using the 
positions and sizes of over 200 000 random close packed 
particles from our experiment, we confirm that our ex- 
perimental system is hyperuniform and incompressible. 
Our results are consistent with prior work [l0[ and the 
data can be used with other algorithms for quantifying 
hypcruniformity in polydispcrsc samples |ll| . We also 
sec a relationship between local polydispersity and lo- 
cal order, confirming that locally a higher polydispersity 
results in less ordered packing. 
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